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ABSTRACT 

This paper analyzes the integrated ultraviolet spectra taken with the International 
Ultraviolet Explorer (lUE) and optical ground-based spectra taken with the William 
Herschel Telescope (WHT) of the giant Hii region NGC 604. These data are comple- 
mented with ultraviolet (WFPC2 through F170W) and Ha (WFPC2 through F656N) 
images retrieved from the Hubble Space Telescope (HST) archive. 

The ultraviolet resonance wind stellar lines, the nebular optical emission lines and 
the higher order terms of the Balmer series and Hel absorption lines detected in the 
spectra of NGC 604 are interpreted using evolutionary models optimized for young 
star forming regions. The evolutionary state and the massive stellar content of the 
region is derived in a self-consistent way. 

The three techniques applied suggest that the central ionizing cluster in NGC 
604 is very young, ~ 3 Myr old, and that the stars in the cluster were formed in 
an instantaneous burst following a Salpeter (a=2.35) or flatter (a=1.5) initial mass 
function (IMF), having stars more massive that 80 M0. The stellar cluster is able 
to provide most of the ionizing photons needed to photoionize the whole nebula, and 
the wind power to form the central shell structure where the cluster core is located. 
The stellar cluster is affected by an extinction similar to the average extinction that 
affects the ionized gas. The estimated number of massive stars in the cluster is also in 
agreement with that derived from previous studies based on the detection of individual 
stars. The results that we present here support the use of these techniques for the 
interpretation of the integrated light of more distant star forming regions. 

Key words: techniques: spectroscopic - galaxies: stellar clusters - ultraviolet: stars - 
ISM: Hii regions - ISM: individual: NGC 604. 



1 INTRODUCTION 

Giant extragalactic H ll regions (RH lis) are amongst the 
brightest objects in galaxies. RHlls have been studied ex- 
tensively because they are the best indicators of the con- 
ditions that lead to massive star formation, and they show 
the cloud properties immediately after the moment when 
stars form. Equally important, they are favourable places 
where to measure the chemical abundances and star forma- 
tion rates in galaxies, and hence can be used to trace the 
chemical and star formation history of the universe, when 
observed in galaxies at different redshift (Madau et al. 1996; 
Kobulnicky & Zaritsky 1999). 

RHlls are characterized to have a size larger than 100 
pc and Ha luminosity brighter than 10'^^ erg s~^ (Kenni- 
cutt 1984). Therefore, the nebula requires an ionizing pho- 
ton luminosity larger than 10^^ s~^; this is provided by a 
stellar cluster that contains more than 100 young massive 



stars. These characteristics are very similar to those of star- 
burst galaxies. However, they are less luminous than pro- 
totypical starbursts, and thus they are referred to as mini- 
starbursts (Walborn 1991). Like starbursts, RHlls show a 
nebular emission-line spectrum at optical wavelengths, and 
an absorption-line spectrum at wavelengths shorter than 
the Balmer jump (Leitherer 1997; Rosa, Joubert & Ben- 
venuti 1984). This spectral morphology reflects the fact that 
RHlls are powered by massive stars. These stars emit pho- 
tons with energies of tens of eV which are absorbed and 
re-emitted in their stellar winds, producing ultraviolet reso- 
nance transitions. However, the stellar wind is optically thin 
to most of the ultraviolet photons, that can travel tens of 
parsec from the star before they are absorbed and photoion- 
ize the surrounding interstellar medium. Subsequently, this 
ionized gas cools down via an emission spectrum. This spec- 
tral dichotomy picture allows to derive the stellar content 
and the evolutionary state of the cluster, through the anal- 
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ysis of the ultraviolet (e.g. Vacca et al. 1995; Leitherer et al. 

1996) , or the optical light (e.g. Cid Fernandes et al. 1992; 
Garcia- Vargas & Di'az 1994; Stasiiiska & Leitherer 1996) us- 
ing evolutionary synthesis and photoionization models. The 
two techniques have been applied to the prototypical star- 
burst nucleus NGC 7714, giving similar results (Gonzalez 
Delgado et al. 1999a). 

However, around the Balmer jump the spectra of star- 
bursts (Gonzalez Delgado et al. 1998) and some Hll regions 
(Terlevich et al. 1996) show the higher order terms of the 
Balmer series and He I lines in absorption, formed in the 
photospheres of massive stars. These stellar lines can be de- 
tected in absorption because the strength of the gaseous 
Balmer lines in emission decreases rapidly with decreasing 
wavelength, whereas the equivalent width of the stellar ab- 
sorption lines is almost constant with wavelength (Gonzalez 
Delgado, Leitherer & Heckman 1999b). Evolutionary syn- 
thesis models that predict the profiles of the higher order 
terms of the Balmer series and He I lines in absorption can 
also be used to estimate the evolutionary state of the stellar 
cluster in Hll regions. 

The second more luminous and brightest H ll region in 
the Local Group of galaxies, after 30 Dor in the LMC, is 
NGC 604 in M 33. Its distance, 840 kpc (Freedman et al. 
1991), allows detailed studies of the individual stars, and 
also to obtain integrated properties. Thus, NGC 604 is an 
excellent laboratory in which to explore questions about the 
effect of star formation, the slope and upper mass limit of 
the IMF, and to test the consistency between the three tech- 
niques described above, that allow to determine the stellar 
content and the evolutionary state of the ionizing stellar 
cluster from the spatially integrated ultraviolet and optical 
spectra of the Hll region. 

NGC 604 has been intensively studied in the past (e.g. 
Peimbert 1970; Israel & van der Kruit 1974; Smith 1975; 
Hawley & Grandi 1977; Kwitter & AUer 1981; Conti & 
Massey 1981; Rosa & D'Odorico 1982; Viallefond & Goss 
1986). A precise determination of the chemical composition 
of the gas was made by Di'az et al. (1987); they give an 
oxygen abundance 12-(-log(0/H)=8.5. VLA observations in- 
dicate that the gas is very tenuous (average electron den- 
sity, rms A^c, between a few and ~ 100 cm~^), and not 
very dusty, with a mean visual extinction Av ~ 0.5 mag 
(Churchwell & Goss 1999). The morphology of the ionized 
gas is very complex, showing many filaments and shell struc- 
tures that are expanding (e.g. Hippelein & Fried 1984; Rosa 
& Solf 1984; Clayton 1988; Sabalisck et al. 1995; Muiioz- 
Tunon et al. 1995; Yang et al. 1995; Medina Tanco et al. 

1997) . This morphology is a consequence of the violent star 
formation activity in NGC 604. HST imaging photometry 
has revealed that the stellar cluster is resolved into ~200 
massive stars in an area of ~ 10000 pc^ (Drissen, Mofi'at & 
Shara 1993; Hunter et al. 1996). Evidences that the stellar 
cluster is evolved come from the detection of Wolf-Rayet 
stars (Conti & Massey 1981; D'Odorico & Rosa 1981; Dris- 
sen, Moffat & Shara 1990, 1993), one supernova remnant 
(D'Odorico et al. 1980), and one candidate to red super- 
giant (Terlevich et al. 1996). The existence of these stars as 
members of the stellar cluster of NGC 604 suggests that the 
age of the region is in the range 3-5 Myr. 

This paper presents spatially integrated ultraviolet and 
optical spectra of NGC 604. The goal is to derive the mas- 



Figure 1. HST+WF'PC2 (F170W) image of NGC 604 at ultra- 
violet wavelengths. The scale is 0.046 arcsec pixel"^, equivalent 
to 0.19 pc pixel" ^. The rectangles mark the positions of the lUE 
apertures. North is up and East to the left. 



Figure 2. F5T-I-WFPC2 (F653N) image in the Ha emission 
line (grayscale). The ultraviolet emission (F170W) is overlayed in 
white contours. The scale is 0.1 arcsec pixel"^. The lines indicate 
the slits of the ground-based optical spectroscopy. Note that the 
core of the cluster is within a nebular hole of 20 pc radius. North 
is up and East to the left. 

sive stellar content and the evolutionary state of the ionizing 
cluster, by means of evolutionary synthesis and photoioniza- 
tion models applied to the ultraviolet resonance wind stel- 
lar lines, to the nebular emission lines and to the higher- 
order terms of the Balmer series and He I lines in absorption. 
The consistency between the results obtained with the three 
techniques, and the similarity with the results obtained from 
studies based on the detection of individual stars, strength- 
ens the reliability and power of these techniques when they 
are applied to determine the stellar content of more dis- 
tant star-forming regions. Section 2 presents the observa- 
tions. Section 3 describes the photoionization models that 
fit the emission-line spectrum of the nebula. In section 4, 
the higher-order terms of the Balmer and He I lines in ab- 
sorption are analyzed. Evolutionary synthesis models of the 
ultraviolet resonance wind stellar lines are in section 5. In 
section 6, the massive stellar content is derived. The sum- 
mary and conclusions are in section 7. 



2 OBSERVATIONS 
2.1 HST images 

NGC 604 has been intensively observed with the WFPC2 
camera onboard the HST. For the purpose of this work, 
we have retrieved HST archive images at the ultraviolet 
(F170W fiher) and optical Ha (filter F656N) wavelengths. 
Two F170W exposures of 350 s each and two F656N ex- 
posures of 1000 s each were combined to produce one final 
ultraviolet (Fig. 1) and one Ha image (Fig. 2). The nebula 
was centred in the PC camera, therefore the ionizing clus- 
ter is sampled at 0.046 arcsec pixel"^, corresponding to 0.19 
pc pixel"'^ (at a distance of 840 kpc, 1 arcsec corresponds 
to 4.1 pc). The ultraviolet emission is spread in the inner 
20x20 arcsec. The stellar density looks much lower than 
that observed in super star clusters. Surface brightness pho- 
tometry indicates that 75 percent of the ultraviolet flux is 
produced by the core of the stellar cluster, that is contained 
within a central nebular hole of diameter ~ 15 arcsec; we 
shall refer to this as the centre of the nebula. The Ha image 
shows many filaments and shell structures that extend out 
to ~ 200 pc from the centre of the nebula. However, the 
brightest parts of the nebula trace almost a ring structure 
at a distance of 40 pc from the centre. 

The total ultraviolet flux per unit wavelength in the 
F170W image is 1.9 x 10"^^ erg s"^ cm"^ A"^ The Ha 
flux is 4.0 X 10^^^ erg s~^ cm~^, in good agreement with 
the value measured by Churchwell & Goss (1999). 
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Table 2. Input parameters to 
CLOUDY. 



12+log(0/H) 
log (He/H) 
log (N/H) 
log (Ne/H) 
log (S/H) 
logQ 

Inner Radius 
Ne 



8.51 

-1.08 

-4.65 

-4.27 

-5.05 

51.54 (s-i) 
20 pc 

30, 100 cm-3 



of the gas to the ionizing source. For a spherical geometry, 
the average U is proportional to {(jP" N^QY^^ ■ 

The chemical abundances are fixed to those derived 
by Diaz et al. (1987). We scale the chemical composi- 
tion of the gas to the value of the oxygen abundance, 
12-|-log(0/H)=8.5, except the abundance of N, Ne, S and 
He, for which we take the values derived by Diaz et al. 
(1987; see Table 2). The electron density derived by Diaz 
et al. (1987) is of the order of 100 cm~'^. Our measurements 
of the [S ll]6717/6731 ratio indicate the low-density limit, 
with values below 100 cm~"^ ^. We have computed constant 
density models at two values of the electron density, 100 
cm~'^ and 30 cm~^. The latter value comes from the filling 
factor and an rms electron density of 10 cm~^, derived from 
the Ha image. Also, an electron density of 30 cm~'^ predicts 
a [S ll]6717/6731 ratio which is in better agreement with the 
observed value. 



3.2 The radiation field 

The radiation field used as input to CLOUDY is the spectral 
energy distribution generated by the evolutionary synthesis 
code developed by Leitherer and colaborators (Leitherer et 
al. 1999). The code includes the new set of stellar evolu- 
tionary models of the Geneva group (Schaller et al. 1992; 
Schaerer et al. 1993a, b; Charbonnel et al. 1993; Meynet et 
al. 1994), and the stellar atmospheres grid compiled by Leje- 
une et al. (1997), supplemented by the expanding spherical 
non-LTE models of Schmutz, Leitherer & Gruenwald (1992). 
The latter stellar atmosphere models are applied to stars 
with very strong stellar winds. The spectral energy distri- 
bution was generated using the Z=0.008 (Z0/2.5) metal- 
licity tracks, assuming that stars and gas have the same 
metallicity. We assume two different star formation scenar- 
ios: instantaneous burst and continuous star formation at 
a constant rate. Different models are also computed with 
different assumptions about the IMF. The slope is Salpeter 
(a= 2.35), flatter (a= 1.5), or steeper (q= 3.3). The upper 
mass limit cut-off is set to 30, 60, 80 or 120 M0. The lower 
mass limit is fixed to 1 M0. However, photoionization mod- 
els are not influenced by the lower limit mass cut-off if it is 
below 10 M0. The spectral energy distribution is normal- 
ized to an ionizing photon luminosity of logQ = 51.54 s~^. 
This value is derived from the total Ha fiux measured in the 
image, after correcting by the mean extinction, Av = 0.5 
mag, derived by Churchwell & Goss (1999). We will check 



S Notice, however, that the atomic parameters for the sulfur have 
changed in the intervening time 



below the compatibility of this normalization with the mass 
of the stellar cluster derived from the ultraviolet continuum 
luminosity. 



3.3 The observational constraints 

The emission lines were measured in our scanned spectrum 
of the nebula. The observed ratios have been corrected by 
the mean extinction, c(H/3)=0.22, derived by Churchwell & 
Goss (1999). This value is in agreement with the redden- 
ing derived from the H7/H/3 ratio measured in the scanned 
spectrum, when this ratio is corrected by the underlying 
stellar absorption. This extinction is lower than the values 
derived by Di'az et al. (1987), c(H/3)=0.3-0.4; however, the 
extinction changes across NGC 604 (Mai'z-Apellaniz 1999). 
The de- reddened emission line ratios to be fitted by the pho- 
toionization models are included in Table ^. This table also 
includes the emission line ratios derived by Diaz et al. (1987) 
in several parts of the nebula, and those from an integrated 
spectrum in Mai'z-Apellaniz (1999). The range of values for 
each of these emission line ratios defines the tolerance that 
we accept for the difference between the observed and the 
predicted values from the models. 

Following Stasiriska & Schaerer (1999), we try to fit the 
strength of the emission lines, and also emission-line ratios 
indicative of the electron temperature ([O lll]4363/5007), 
electron density ([S ll]6717/6732) and ionization structure 
([Oiii]5007/[Oii]3727, [Siii]9069/[S ii]6716+6732). Another 
important observational constraint is the radius of the neb- 
ula. The Ha image shows that the flux extends out to 200 
pc of the center of the region; however, more than 90, 75 and 
50 percent of the total flux is within the inner 140 pc, 100 
pc and 60 pc radii, respectively (Fig. 3). Our determination 
of the core radius of the region is 80 pc. 



3.4 Model results 

We start by fixing the filling factor and the electron den- 
sity by means of fitting the line ratios [S n]6717+6731/H/? 
and [S ll]6717/6731, and the radius of the nebula. 
[S ll]6717+6731/H/9 is a good calibrator of the ionization 
parameter, because this ratio does not depend much on 
the IMF assumptions, on the star formation law or on the 
evolutionary state of the cluster (Gonzalez Delgado et al. 
1999a). Instead, this ratio depends on the filling factor, and 
for fixed geometry, ionizing photon luminosity and electron 
density, it depends on the ionization parameteiPj]. The ob- 
served ratio, [Sll]6717+6731/H/3 ~ 0.4, indicates a filling 
factor <;!) ~ 0.1 (Fig. 6). Therefore, a large fraction of the 
volume of the region is filled with ionized gas; however, this 
gas is very tenuous because the rms electron density mea- 
sured in the Ha image is ~ 10 cm~"^, and the electron density 
in the S"*" zone, derived from the definition of filling factor 
{(j, = {N4rms)/Nef), is ~ 30 cm . The observed ra- 
tio [S ll]6717/6731=1.4, indicates that the electron density is 
< 100 cm^^*. Models with a filling factor of 0.1 and Na = 30 

^ Nonetheless, we should always be cautious when using sulfur 
dependent diagnostics, given the uncertainties about the atomic 
parameters and its not so well known ionization; also this calibra- 
tion relies on the assumption of an ionization bounded nebula. 
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Figure 7. Predicted emission-line ratios, rj parameter and the electron temperature as a function of the cluster age, for an instantaneous 
burst (filled symbols) and continuous star formation (open symbols). The models assume that the mass of the cluster is distributed 
following a Salpeter IMF between 1 and 120 M0. The filling factor is 0.1 and the electron density 30 cm~'^. The observed values are 
represented by an horizontal line. 



4.1 Description of the models 

Gonzalez Delgado et al. (1999b) have computed evolutionary 
stellar population synthesis models that predict the photo- 
spheric absorption H Balmer and He I lines, between 3700 
and 5000 A, for a single-metallicity stellar population. The 
models, which are optimized for galaxies with active star 
formation, synthesize the profiles of the H Balmer series 
(H/3, H7, H5, H8, H9, HIO, Hll, H12 and H13) and the 



He I absorption lines (He I A4922, He I A4471, He I A4388, 
He I A4144, Hel A4121, Hel A4026, Hel A4009 and Hel 
A3819), with a spectral sampling of 0.3 A pixel"^, for a 
burst and for continuous star formation at a constant rate. 
They use a stellar library that includes NLTE absorption 
profiles for stars hotter than 25000 K, and LTE profiles for 
lower temperatures. The temperature and gravity coverage 
is 4000 < Teff < 50000 K, and 0.0 < log g < 5.0, respectively 
(Gonzalez Delgado & Leitherer 1999). The models assume 
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Figure 8. As Figure 7, but for an instantaneous burst with a Salpeter IMF upper mass limit cut-off 120 M0 (filled circles), 80 M0 
(open circles), 60 M0 (squares), and 30 M0 (open triangles). 



that stars evolve from the main sequence following the evo- 
lutionary tracks of the Geneva group (Schaller et al. 1992; 
Schaerer et al. 1993a,b; Charbonnel et al. 1993; Meynet et al. 
1994). The strength of the Balmer and He I lines is sensitive 
to the age after the first 3 Myr of evolution, and sensitive to 
the IMF if the age is younger than 3-4 Myr. Models assume 
that stars have a metallicity Z — 0.008, and make different 
assumptions about the IMF slope (a=1.5, 2.35 and 3.0), and 
upper mass limit cut-off (Afup = 80 M0 and 30 M0). For 
the purpose of this paper, the range of ages computed spans 
from to 10 Myr. 



4.2 Model Results 

We use the equivalent width of He I A4388, A4026, and 
A3819, and the wing absorption profiles of H/3, H7, H5, H8, 
H9 and HIO to constrain the age, IMF and star formation 
law. He I A4388, A4026, and A3819 are detected in absorp- 
tion, because the corresponding nebular emission lines have 
equivalent widths that are at least a factor 3 weaker than 
the nebular emission in He I A4471 (Gonzalez Delgado et al. 
1999b). Note that He I A4471 is only partially filled with 
the nebular emission; therefore, the nebular emission of the 
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lines He I A4388, He I A4026 and He I A3819 is ~0, and the 
equivalent widths of these absorption features represent well 
the strength of the stellar continuum radiation of the clus- 
ter. The equivalent widths of these lines, He I A4388, He I 
A4026 and He I A3819, measured in the spectrum of Fig. 
5 are 0.32, 0.55 and 0.26 A, respectively. Their strength is 
compatible with continuous star formation and with burst 
models 3 Myr old (Fig. 10). However, He I lines can con- 
strain the IMF. Models with an IMF steeper than Salpeter 
or with Mup = 30 M© predict He I lines stronger than ob- 



served. The strengths of the lines are compatible with a 3 
to 4 Myr old burst formed following a Salpeter or slightly 
flatter IMF (a= 1.5, Fig. 10). The absorption wings of the 
Balmer lines are also compatible with burst models 3 Myr 
old, with Salpeter or slightly flatter IMF (a= 1.5, Fig. 11). 
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burst galaxie^ Most of these lines are formed in the stellar 
winds of the massive stars that belong to the starburst. The 
profile of these lines reflects the stellar massive content in 
the starburst; therefore, they depend on the IMF and star 
formation law (Leitherer, Robert & Heckman 1995). In this 
section, we constrain the evolutionary state and the IMF of 
NGC 604 by means of fitting the profiles of the ultraviolet 
stellar lines Siiv and Civ. 

5.1 Description of the models 

Evolutionary stellar population models have been computed 
with the code Starburst 99 (Leitherer et al. 1999). The code 
uses a stellar library of lUE ultraviolet spectra of O, B and 
Wolf-Rayet stars (Robert, Leitherer & Heckman 1993). The 
spectral resolution of the O and Wolf-Rayet stars observed 
is 0.75 A. These stars are located in the solar neighborhood. 
However, the evolutionary models computed here assume 
that the stars evolve from the main sequence following the 
evolutionary track at .Z = 0.008. Models are computed for 
instantaneous bursts between and 10 Myr, and for continu- 
ous star formation lasting 10 Myr, and different assumptions 
about the slope (a=2.35, 3.0 and 1.5) and upper mass limit 
cut-off (Mup =120, 80, 60 and 40 M0) of the IMF. The 
model spectra are smoothed to the lUE spectral resolution 
of the observations of NGC 604, which is 6 A. 

5.2 Model results 

The strongest wind stellar features in the spectra of NGC 
604 are Nv A1240, Siiv A1400, Civ A1550, Hen A1640, and 
Niv A1720. These features are present in all the spectra of 
NGC 604; however, for this analysis we use the average spec- 
trum of those at PA=149°, 136° and 159°, and the spectrum 
at PA=32°, because they have the best signal to noise ratio. 
The evolutionary state, the star formation law and the IMF 
in NGC 604 is constrained with the profiles of these lines. 

Hell shows a broad emission profile in Wolf-Rayet and 
03-05 supergiant stars; N iv shows a PCygni profile in 
Wolf-Rayet stars; N v and C iv show strong PCygni profiles 
in all O stars, and Si iv only in O supergiants. Thus, the pro- 
file of these lines in the integrated spectrum of a starburst 
depends strongly on the stellar content and age of the stellar 
cluster. In particular, Siiv shows a strong PCygni profile if 
the cluster formed in an instantaneous burst and its age is 
between 3 and 5 Myr, because within this age interval the 
ultraviolet light is dominated by O blue supergiants. On the 
other hand, the profile of C iv depends strongly on the IMF. 
It shows a strong PCygni profile if stars more massive that 
60 M0 are formed in the cluster, and if the slope of the IMF 
is flatter than Qf=3.0. 

The spectrum of NGC 604 indicates that its stellar pop- 
ulation must be dominated by massive, young O stars; thus, 
the stellar cluster must be young and of short duration. Con- 
tinuous star formation models can be ruled out because they 
show Siiv weaker than observed (Fig. 12). In these models 
the line is diluted because the fraction of O supergiants with 

II One important exception is knot A in the RHn NGC 2363 
(Drissen et al. 2000), so young that the effect of winds is not yet 
reflected in its uv spectrum. 



respect to the total number of O stars is lower in the contin- 
uous star formation than in the instantaneous burst models. 
The strength of Siiv in NGC 604 indicates that it is an in- 
stantaneous burst. The age of the burst has to be between 3 
to 5 Myr, because bursts younger than 3 Myr or older than 
5 Myr have very few O supergiants, and thus very weak Si iv 
(Fig. 13) . On the other hand, from the strength of C iv we 
can exclude instantaneous bursts with Mup < 60 M0 (Fig. 
14), and IMF steeper than a=3.0 (Fig. 15). Therefore, we 
conclude that continuous star formation and instantaneous 
burst models with few very massive stars can be ruled out. 
The wind lines are compatible with an instantaneous burst, 
formed following a Salpeter or slightly flatter IMF, with up- 
per mass limit cut-off higher than 60 M0. A 3 Myr instan- 
taneous burst with Salpeter and Mup > 80 M0 fits well the 
profile of the wind stellar lines. 



6 THE MASSIVE STELLAR CONTENT OF 
NGC 604 

We derive the massive stellar population in NGC 604 from 
the ultraviolet continuum luminosity, comparing the obser- 
vations with the predictions of synthetic models. The num- 
ber of Wolf-Rayet stars are derived from the luminosity of 
the He ll A1640 line. The stellar content derived in this way is 
compared with that derived from HST photometry of stars 
in NGC 604. However, we need first to estimate the extinc- 
tion. 

6.1 Extinction estimates 

Leitherer & Heckman (1995) have shown that the ultraviolet 
continuum arising from a young starburst (age < 10 Myr) 
has a spectral index, /? (Fa oc A''), which is independent 
of the IMF, star formation law and metallicity. Thus, any 
deviation from the predicted value, ~ 2.5 for a 3 Myr in- 
stantaneous burst, could be attributed to reddening. 

The ultraviolet continuum flux distribution of the lUE 
spectra of NGC 604 shows evidence of reddening, because af- 
ter correcting by Galactic extinction, E{B — V) = 0.03 (Mc- 
Clure & Racine 1969) , the spectra are flatter than the spec- 
tral energy distribution predicted by the evolutionary syn- 
thesis models. Massey & Hutchings (1983), analysing lUE 
spectra of H ll regions in M33, conclude that the M33 extinc- 
tion curve is significantly different from that of the Galaxy, 
because the spectra of the H ll regions show very weak 2200 
A interstellar absorption dip, resembling those of the LMC 
and SMC. Therefore, we use the LMC curve to match the 
observed spectra with the spectral energy distribution of a 
3 Myr old instantaneous burst. The E{B — V) derived is 
0.1 for the spectrum at PA=32°, and 0.12 for the average 
spectrum of those at PA=136°, 159° and 149°. These values 
are in agreement with those derived by Massey & Hutchings 
(1983). However, the extinction can be 0.03 higher than the 
values above if it is derived comparing the observed spec- 
tra with the spectral energy distribution predicted includ- 
ing only the stellar contribution. This is due to the nebular 
contribution being dominated by the two-photon continuum 
emission that peaks at ~1500 A. If the nebular continuum 
is not included, the spectral energy distribution is steeper 
and a higher extinction is required to match the observed 
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with that derived from photometric studies based on the 
detection of individual stars. 
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